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We present a comprehensive study of the temperature and doping dependence of the 500 meV 
peak observed at q = (7r,0) in resonant inelastic x-ray scattering (RIXS) experiments on La2Cu04. 
The intensity of this peak persists above the Neel temperature (Tjv=320 K), but decreases gradually 
with increasing temperature, reaching zero at around T=500 K. The peak energy decreases with 
temperature in close quantitative accord with the behavior of the two-magnon Big Raman peak 
in La2Cu04, and with suitable rescaling, agrees with the Raman peak shifts in EuBa2Cu306 and 
K2NiF4. The overall dispersion of this excitation in the Brillouin zone is found to be in agreement 
with theoretical calculations for a two-magnon excitation. Upon doping, the peak intensity decreases 
analogous to the Raman mode intensity and appears to track the doping dependence of the spin 
correlation length. Taken together, these observations strongly suggest that the 500 meV mode is 
magnetic in character and is likely a two-magnon excitation. 

PACS numbers; 78.70.Ck, 78.30.-j, 71.10.Fd, 75.10.Pq 



I. INTRODUCTION 

Since the discovery of high temperature supercon- 
ductivity, the cuprate materials have provided a fertile 
ground for studying the physics of electron correlations 
and quantum magnetism. Various spectroscopies have 
been instrumental in revealing the electronic struc- 
ture of these materials, ^"'^ mostly focusing on the low 
energy region below 100 meV, which is of order of 
the pseudogapi^ However, one aspect that separates 
cuprates from many other condensed matter systems 
is their intrinsically large energy scale - a result of the 
strong hybridization between the in-plane copper and 
oxygen orbitals4 This fact has driven recent interest 
in spectroscopic studies of the high energy region. 
Examples include neutron scattering studies of the 
universal spin excitation spectrum;^^— optical studies 
of the well-known mid-infrared featuresji^r— and angle 
resolved photoemission investigations of high energy 
kinks. Complementing these methods is a relatively 
new technique, resonant inelastic x-ray scattering 
(RIXS), which provides momentum resolved information 
on electronic excitations over a wide energy rangeJ^^— 

A recent Cu i^-edge RIXS study of insulating and 
doped La2-xSrxCu04 and Nd2Cu04 found a new exci- 
tation at 500 meV at a momentum transfer q of (tt, 0).'^^ 
This mode was found to soften and lose its spectral 
weight away from the (tt, 0) position, and to weaken 
in intensity upon doping holes into the copper-oxygen 
plane. The energy scale and doping dependence of this 



mode point towards a magnetic origin, and it was argued 
that the excitation is related to the two-magnon which 
has been observed previously with Raman scattering'^®^— 
and optical absorption4^"— While a previous RIXS 
study at the oxygen Is resonance in Sr2Cu02Cl2 showed 
a weak feature at 0.5 eV that was associated with a 
two-magnon,'*'' this comprehensive Cu i^-edge study^ 
opened an exciting new avenue of studying magnetic 
excitations with the RIXS technique. Subsequently, exci- 
tations with a two-magnon-like dispersion were observed 
in both cuprates and NiO with high-resolution L-edge 
RIXS^-^ and M-edge RIXS.^* Although excitations 
in La2Cu04 at the zone boundary were not presented 
in this soft x-ray study, the extrapolated dispersion 
suggests strong similarity with the hard x-ray RIXS 
result, giving further credence to the two-magnon in- 
terpretation. Nevertheless, other possibilities, such as a 
d-d crystal field excitatiou )^^'^^ could not be completely 
ruled out as possible explanations for the 500 meV mode. 

In order to further understand this mode, we have 
carried out comprehensive RIXS studies of the temper- 
ature, doping and momentum dependence of the 500 
meV peak. We observe that the RIXS peak becomes sig- 
nificantly weaker in intensity and shifts to lower energy 
as the temperature is raised above room temperature, 
and the peak appears to vanish at around 500 K. This 
behavior is very similar to the temperature dependence 
of two-magnon Raman scattering in YBa2Cu306 and 
EuBa2Cu306 studied by Knoll et al<^ The momentum 
dependence of the peak energy and spectral weight in the 
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insulator is found to be in semi-quantitative agreement 
when compared with recent theoretical calculations of 
the dispersion of two-magnon excitations.—"— Finally, 
hole doping causes the peak to shift to lower energy and 
to lose intensity in a manner very similar to the behavior 
of the two-magnon Raman peaki^ Taken together, the 
above observations strongly suggest that the observed 
500 meV peak corresponds to a two-magnon excitation. 

Despite good agreement in temperature and doping 
dependence, there remains an important difference 
between the two-magnon Raman and the RIXS data. 
As pointed out in Ref. HI, a surprising result was that 
the 500 mcV RIXS peak energy is significantly higher 
than the ~400 meV low-temperature 2-magnon Raman 
peak.— 1^ Theory predicts that the 2-magnon Raman 
peak energy should be renormalized to around 2.75J by 
strong magnon-magnon interactionsj^^'^'^ where J is 
the effective magnetic superexchange energy including 
the Oguchi correction. Such an estimate, for which J 
145 meV, would place the 500 meV RIXS peak energy 
slightly below 3.5J. Since 3.5 J corresponds to the peak 
in non- interacting density of states for q=(7r 0), this 
implies that magnon-magnon interaction is much less 
than the q=(0 0) case probed by Raman scattering)^ 
One possible explanation given in Ref. [s^ was that 
the two spin flips are created on adjacent copper-oxide 
planes. Here, we suggest another possibility — namely, 
that the interaction between magnons in the same plane 
is weakened for total momentum q — (tt, 0) (as compared 
to q = (0,0) which is explored in Raman scattering). 

The paper is organized as follows. Section |ll] out- 
lines the experimental method. Section IIIII presents 
the results for the incident energy and momentum 
dependence, the temperature dependence of the mid-IR 
peak and the effect of doping with holes. In Section 
IIV A[ the temperature dependence is discussed in 
more detail, compared to other materials, and qual- 
itatively described with a simple model. In section 
IIVB| we discuss a possible reason for the apparently 
non-interacting value of the observed peak position. 
Section |V] contains the summary and our conclusions. 
The various background subtraction methods used in 
analyzing the data are described in detail in Appendix A. 



II. EXPERIMENTAL DETAILS 

The same La2Cu04 sample used in Ref. was used 
for all of the La2Cu04 measurements presented here. 
We also studied underdoped Lai.93Sro.o7Cu04 to further 
examine the doping dependence of the 500 meV feature. 
This crystal was grown at the University of Toronto 
by the traveling solvent floating-zone method, and cut 
along the ac plane. After annealing in flowing oxygen 
for 30 hours at 900 °C, the superconducting transition 
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FIG. 1: The scattering configurations used when scattering 
in the (a) vertical (sector 9ID experiments) and (b) horizon- 
tal plane (at sector 30ID and BLllXU). The arrows indicate 
the wavevectors of the incident and scattered beams. The 
sample is shown in both cases along the y-axis of the figure 
(note that the vertical direction in real space corresponds to 
perpendicular to the incident beam, and parallel to the plane 
of the page in (a), while it is perpendicular to the plane of 
the page in (b)) The crystallographic axes are also indicated. 
In (a), the incident beam polarization is directed out of the 
page, parallel to the c-axis. In (b), the incident polarization 
is in the scattering plane. The scattering angle 2Q is set to 
nearly 90°, which suppresses the in-plane polarization of the 
scattered beam. The angle & is about 20°, so as to have a 
large c-axis component of the incident polarization. 



temperature of this sample was 14 K, as determined by 
the onset temperature of the diamagnetic signal. The 
RIXS experiments were carried out at the XOR-IXS 
30ID and 9ID beamlines of the Advanced Photon Source, 
and at the BLllXU beamline at SPring-8, Japan. The 
energy calibration was identical for all of the beamlines, 
with the absorption iiT-edge of a copper foil set to 
8980.5 (NP The measurements were all carried out by 
flxing the incident energy, Ei , and scanning the scattered 
photon energy, Ej. Data are plotted as a function of 
energy loss uj = Ei — Ef. 

The RIXS measurements at the 30ID beamline were 
performed using the MERIX spectrometer. X-rays 
impinging upon the sample were monochromatized to a 

bandwidth of 72 meV, using a four-bounce (-1 h) 

monochromator with asymmetrically cut Si(400) 
crystals. The beam size on the crystal was reduced 
to 45(H)x20(V) ^m^ by two segmented piezoelectric 
bimorph mirrors, in a Kirkpatrick-Baez conflguration 
(ACCEL Instruments). The photon flux on the sample 
was 1.2 X 10^^ ph/s. A Ge (337) spherical diced 
analyzer and position sensitive microstrip detector (with 
channels spaced 125 /xm) were placed on a Rowland 
circle of 1 m radius, in order to select Ef. Using a 
microstrip detector rather than a point-detector allowed 
measurement of multiple energy channels simultaneously 
and reduced geometrical broadening of the spectral 
resolution functioni^ As the center channel energy 
is varied during the energy scan, the same energies 
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are measured multiple times and binned at the end 
of the scan. The overall energy resolution for the 
measurement, with the sample in place, was 120 nieV 
(FWHM) as determined by the width of the elastic peak. 

The scattering plane was horizontal, with the in- 
cident polarization parallel to the scattering plane 
(TT-polarization) and the scattering angle was set to 
nearly 90 degrees, thereby suppressing the background 
from the elastic peaki^i^li^dliS^ In order for the 
incident beam polarization to have a large c-axis com- 
ponent (as in Ref. [ssl ) while maintaining a ^90 degree 
scattering angle, the measurements were carried out at 
the Q-positions of (3.5 6), (3.5 0.25 6), (3.5,0.5,6) 
and (3 8.5) for reduced momentum q of (tt 0), (tt 
7r/2), (tt tt) and (0 0), respectively. Here the wavevector 
change of the scattered beam Q = q + G, where G is 
the reciprocal lattice vector closest to Q. We use the 
tetragonal notation, for which the lattice parameter is 
~3.78 A along the Cu-O-Cu direction. A diagram of the 
scattering geometries is shown in Fig. [TJ Momentum 
and incident energy dependence measurements were 
performed at room temperature, while a closed-cycle He 
refrigerator was used for the temperature-dependence 
study, which was carried out in a temperature range of 
35 K to 500 K. 

The experimental condition at the 9ID beamline, 
scattering in the vertical plane, was identical to that of 
Ref.m with an overah resolution of 130 meV (FWHM). 
All of the measurements at 9ID were carried out near 
the Q-position of (2.5 0), which corresponds to the 
reduced momentum q of (tt, 0). Q=(3 0) was used 
for zone center momentum transfer, and (2.5 0.5 0) was 
used for the q=(7r, tt) equivalent position. The sample 
was mounted on an Al sample holder and measure- 
ments done in the temperature range from 45 K to 300 K. 

An additional temperature-dependence experiment 
was carried out at the BLllXU beamline of SPring-8 
at the same Q. This instrument had lower resolution 
(440 meV) and used the same horizontal scattering 
configuration described above. Although the lower reso- 
lution in this latter setup prevented us from extracting 
detailed energy parameters, intensity information could 
be obtained from this experiment, and this is reported 
in Sec. IIIIBI The background subtraction methods used 
in the various measurements are explained in detail in 
Appendix A. 
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FIG. 2: (Color Online) (a) Plot of RIXS intensity vs. en- 
ergy loss, measured at various incident energies at the 30ID 
beamline (horizontal scattering geometry), for q=(7r, 0). The 
curves are offset along the y-axis as indicated by the thick grey 
lines on the right. The solid lines are guides to the eye. (b) 
Background-subtracted intensity is plotted as a function of 
energy loss, measured at the 30ID beamline at room temper- 
ature for incident energy of 8994 eV, for various momentum 
transfers, (c) same as (b) but in vertical scattering geometry, 
from Ref. |33 (T=45K). The solids lines are results of fits to 
Lorentzian lineshapes. 



the absorption K-edge in Cu. The solid lines are guides 
to the eye. The incident energy was thereafter set to 
8994 eV, around the peak of the resonance. The 500 
meV excitation resonates at the same photon energy as 
the higher energy 3.5 eV excitation at this q<^ The mo- 
mentum dependence of the energy-loss spectra measured 
at this incident energy is presented in Fig. [5fb)-(c). 
Fig. m^b) shows the room-temperature data taken at 
the 30ID beamline in the horizontal scattering mode, at 
various q. The momentum dep endence is in agreement 
with that observed in Ref. l35l which was measured in 
a vertical scattering geometry. These latter data are 
shown in Fig. [2jc) for comparison. In both data sets, 
the 500 meV peak shows the strongest intensity at 
q=(7r, 0) and disperses to lower energies away from that 
maximum. The peak disappears (or nearly so) at the 
other two high-symmetry points (0,0) and (7r,7r). 



III. EXPERIMENTAL RESULTS 
A. Incident Energy and Momentum Dependence 

Fig. (SJa) shows several RIXS spectra for q=(7r,0) 
obtained with different incident photon energies near 



The main difference between the results from the 
horizontal and vertical geometry measurements is that 
there appears to be some intensity at q=(7r,7r) in the 
horizontal scattering data set. Fig. [H^b), whereas this 
was negligible in the vertical scattering experiment. The 
discrepancy might be due to a non-resonant scattering 
contribution which would remain after the energy-gain 
side background subtraction procedure, as discussed 
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in Appendix A. Polarization could be another possible 
explanation, since in the horizontal scattering case, the 
incoming polarization has a non-zero a-axis component, 
and the outgoing photon has a reduced c-axis compo- 
nent, when compared to the vertical scattering geometry. 
However, a more likely reason for this discrepancy is the 
momentum resolution effect. The momentum resolution 
at q=(7r,7r) in the horizontal scattering experiment 
extended ^--^20% towards (tt 0), while the q-resolution 
in vertical scattering extended only ~6% towards (tt 
0) (the resolution is shown as the horizontal error bars 
in Fig. [3]). The lower momentum resolution means 
that the data obtained at nominal q=(7r,7r) contains 
significant contributions from the neighboring q, giving 
rise to non-zero intensity. When the intensity is a 
rapidly varying function of q, such as may be the case 
near q=(7r,7r),'^^ full resolution deconvolution would 
be required to describe our intensity observations 
quantitatively. Because of the higher q-resolution, and 
a background-subtraction which removes non-resonant 
scattering contributions, we are more confident in the 
vertical scattering result at q=(7r,7r). 

The measured spectra were each fit to a Lorentzian 
plus constant background, shown by the solid lines in 
Fig. [Hb)-(c). The resultant peak positions are plotted 
in Fig. EJa). In Fig. IS^b), we plot the spectral weight, 
calculated by multiplying the Lorentzian peak width 
and intensity parameters together. The spectral weights 
in Fig. mb) are scaled such that the spectral weight at 
q=(7r, 0) is set to unity. Since there is no observable 
Lorentzian peak at q=(0 0), the spectral weight there 
was set to zero (and no error-bars were assigned to those 
zero points). 
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FIG. 3: (Color Online) (a) Dispersion relation of the two- 
magnon excitation obtained from two separate measurements 
as described in the figure legend, and the lines represent the- 
oretical dispersions calculated by Ament et. al.-^^, Donkov 
et al."^ (Aig - with and without final-state magnon-magnon 
interaction included in the calculation, as shown by the dash- 
dotted and vertical dashed lines respectively), Vernay et al.— 
(Aig), and Nagao et al.^^ The energy scale is normalized by 
the magnetic exchange constant J («145 meV). Here we sup- 
pose that 3.5J corresponds to 500 meV. The horizontal er- 
ror bars span the q-resolution. (b) Spectral weight of the 
background-subtracted data, defined as the product of the 
Lorentzian peak intensity and linewidths; this does not in- 
clude constant background, so that the spectral weight at 
zone center is zero. The weights are normalized by the weight 
at q=(7r, 0). Same legend as in (a), two available theoretical 
calculations by Ament et al. and Donkov et al. are plotted. 



Motivated by the original observation of the 500 meV 
mode, there have been a number of recent theoretical 
calculations^"— for the dispersion of the two-magnon 
peak in a Heisenberg antiferromagnet. Nagao and 
Igarashi^ and the van den Brink group ^'^'^^i^^ explicitly 
included a core-hole interaction in their cross-section 
in order to account for the effects of the intermediate 
state in the RIXS process. Donkov et al;^ and Vernay 
et ali^ calculated a generalized momentum-dependent 
Raman response for both Aig and Big scattering 
configuratations, without including core-hole effects. 
A comparison of the theoretical calculations and the 
experimental data is provided in Fig. |31 Since the 
lineshapes of the calculated spectra are typically not 
Lorentzian, we chose to compare the first moment of the 
computed spectra, which is likely to be a more robust 
feature of the theoretical calculations. For clarity, we 
include only the latest calculation from the van den 
Brink group and only the Aig geometry results from 
Refs. [Sa and[56l Although we cannot make a direct link 
between the scattering geometries of RIXS and Raman 
spectroscopies due to the fundamental differences in 
the processes, the Big mode has high intensity and 



energy at q=0, which we do not observe in any our 
RIXS data, so we compare only to the Aig calculations. 
The theoretical energies, originally given in units of 
J, are multiplied by the value of J to compare with 
the experiment. The value of J used here is 143 meV, 
so as to give 3.5J=500 meV. Given the considerable 
differences between the Raman scattering and RIXS 
processes, one should compare these results with caution. 

Fig. [2Ka) shows reasonable agreement between the 
various theoretical dispersion calculations and the 
experiment, especially the RlXS-specific calculations. 
Discrepancies include the fact that when interactions are 
included in the calculation of Donkov and Chubukov,— 
the predicted energies are systematically lower than the 
observed values. As pointed out by these authors, it 
is possible that the random-phase-approximation-type 
analysis overestimated the effect of the interactions. 
The calculated spectral weights are plotted for two of 
the theoretical calculations^^^ in Fig. [3l[b). A main 
point of agreement between these calculations and 
the experiment is the suppression of intensity at the 
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zone center. However, for q=(7r tt) there is a large 
difference between these calculations. This also happens 
to be the q position where the horizontal and vertical 
scattering measurements disagree, but the experimental 
considerations discussed above would support the zero 
intensity at q=(7r tt) result. 

Besides the suppression of spectral weight near zone 
center and relative flattness near q=(7r 0), there is no 
obvious overlap in spectral weight between the theories 
and experiment. This is not very surprising considering 
the many factors that may influence spectral weight as 
a function of momentum, limiting the spectral weight's 
use as a detailed comparison with theory. In contrast 
with this uncertainty, the energy dispersion along (tt 
0)-(7r tt), as seen by the (tt 0) point and the new (tt 7r/2) 
point in Fig. [3][a), is quite consistent with the recent 
RIXS calculations, as are the energies at the other 
q-positions. On the basis of this result it would seem 
that the energy dispersion calculations incorporating the 
core-hole potential^i^i^ agree closely with our data, 
while the pure Raman calculations^!^ agree somewhat 
less consistently. 



B. Temperature dependence 

The 500 meV peak at (tt, 0) was measured at various 
temperatures. The spectra obtained at 30ID are shown 
in Fig. Hl^a) , with a clear 500 meV peak at low tempera- 
tures. As the temperature is increased, the peak becomes 
weaker, and at T=500 K there appears to be little trace 
of the peak left. The data at 500 K can be well described 
as the tail of a Lorentzian centered at zero energy, 
and therefore we use this as our background. At each 
temperature, the Lorentzian fit to the T=500 K data 
was scaled to match the energy loss cj=0 intensity, and 
used as a background (shown as solid lines in Fig.|Ha)). 
Note that for T=35 K, the background calculated in 
this way is below the experimentally observed signal, 
which might imply considerable constant intensity in 
addition to the main peak, as seen in the figure. We do 
not understand the source of this additional intensity 
at present, and here we only focus on the well-defined 
peak. The background-subtracted data are shown in 
Fig. |4] (b)-(d). The solid lines are results of fits to a 
simple Lorentzian function, with an additional constant 
background for the T=35 K data. Given the error 
bars on the data, and that the peak width is not much 
larger than the instrumental resolution, we could not 
reliably extract the temperature dependence of the peak 
widths from the data. The fits were therefore carried 
out simultaneously for all curves keeping the peak width 
the same for all temperatures. The resulting width was 
found to be 166±10 meV. A rapid suppression of the 
intensity as temperature is increased is clearly observed. 
Measurements of the temperature dependence were also 
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FIG. 4: (Color Online) (a) Temperature dependence of (tt, 
0) RIXS spectra (black squares) of La2Cu04 obtained at 
30ID. The background (thin line) was determined by scaling 
a Lorentzian fit of the 500 K data by the elastic line intensity 
at each temperature, as discussed in the text. The spectra 
are displaced along the y-axis for clarity. The thick grey ticks 
on the right side of the graph indicate zero levels, (b)-(d) 
Background-subtracted data and fits for T=400 K, 250 K, 
and 35 K, respectively, (e)-(f) The spectra measured at 9ID 
for temperatures of 250 K, 150 K and 45 K, respectively. The 
background for the 9ID data was obtained from olT-resonance 
data as described in Section [III 



carried out at 9ID and at BLllXU. In Fig. |4Ke)-(g), 
we plot the background-subtracted data obtained at 
9ID. Note that in this case, the off-resonance scans were 
used as background scans as described in Appendix 
A. The solid lines are curves of Lorentzian fits, with 
widths of 153±10 meV which is similar to the 30ID result. 

Figure [5] shows the parameters obtained from the 
Lorentzian fitting of the data as a function of temper- 
ature. Note that the intensity at 500 K comes from 
the residual intensity between 0.4 and 0.5 eV, after sub- 
tracting the fit curve, but is zero within error bar. Also 
shown is the intensity data obtained at BLllXU. The 
Lorentzian linewidth of ~160 meV (FWHM) is close to 
being resolution-limited in all of the scans, and thus it 
was not possible for us to detect any temperature depen- 
dence to the line broadening, although we cannot rule 
it out. The intensity data from the 9ID and BLllXU 
experiments, shown in Fig. [IKa), are scaled to approxi- 
mately match the 30ID intensities at the temperatures 
they share in common. The scatter in the intensity data 
is fairly large but the overall trend clearly shows the in- 
tensity decreasing, reaching nearly zero at 500 K. The 
peak position, plotted in Fig. [HJb), shifts to lower fre- 
quency as temperature increases. The solid line comes 
from fitting to a phenomenological equation, presuming 
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neutron scattering study of the 1-magnon dispersion in 
La2Cu04 has reported that the zone boundary magnon 
energy at T=300 K is about 4% lower than its value 
at T=10 which is also consistent with the current 
result if we assume that the two-magnon energy is 
proportional to the one-magnon energy. Thus, for room 
temperature and below, the shift in the 500 meV peak 
agrees quantitatively with the two-magnon Raman result 
and also scales with the zone boundary 1-magnon energy. 

For higher temperatures, similar temperature depen- 
dence of the excitation energy has been observed for the 
two-magnon Raman spectra of EuBa2Cu306'^ The ex- 
citation energy of EuBa2Cu306 two-magnons obtained 
from the maxima of the fitted curves in the Fig. 1 of Ref. 
[5^ are plotted in Fig. ^h), and compared with that of 
La2Cu04. There is similarity in the temperature depen- 
dence of the two cuprates, although the peak position 
decreases more rapidly with temperature in La2Cu04. 
The temperature dependence will be discussed in more 
detail in Sect. HvXl 



C. Doping Dependence 



FIG. 5: (Color Online) Temperature dependence of the 
Lorentzian fit parameters for tlie RIXS peaks measured at 
q=(7r, 0). Shown are the fit parameters (a) intensity and (b) 
excitation energy. Different symbols correspond to data taken 
with different experimental conditions: horizontal scattering 
at 30ID (Open circles and filled square. The filled square cor- 
responds to a seperate beamtime); Vertical scattering at 9ID 
(filled diamonds); and horizontal scattering at BLllXU (filled 
triangle). The solid line is a result of the fit to Eq. lU, with 
fit parameters So=0.50 eV, 7=1.7, and n=0.1 eV. The dot- 
ted vertical line indicates the Neel temperature for La2Cu04. 
The dashed curve, which is referred to the right-hand scale, is 
the temperature dependence of the two-magnon Raman peak 
position of EuBa2Cu306 obtained from the peaks of the fit- 
ted curves in Ref. [53l . On the same scale, the crosses are the 
two-magnon Raman peak position for La2Cu04 as measured 
in Ref. [371 . 

that a two-magnon mode would interact with ambient 
Bose particles (magnons or phonons), 

EiT)=Eo-{l-j-nin,T)) (1) 

where Eq is the peak energy at T=0 K. 7 and il are 
fitting parameters, and rt(rj,T)= ^^^^^^^^y^_-^ is the 
Bose factor of characteristic energy il. The peak position 
at room temperature is about 5% smaller than at 35 K. 
The magnitude of the peak shift lies on the same curve 
as the shift of the two-magnon Raman scattering peak 
of La2Cu04 observed by Sugai, et. alM, who report a 
3% downshift for T=273 K, shown in Fig. Efb). Note 
that this is also consistent with the Raman peak shifts 
more recently measured by Sugai and Hayamizuj^ (not 
shown in Fig.[SJb) to avoid over-crowding the graph). A 



In addition to the samples measured in Ref. Issl . the 
RIXS spectrum of a La2_xSrxCu04 sample of a;=0.07 
was also measured at 9ID using the vertical scattering 
geometry (the spectrum is shown in Appendix A). 
Combining data from Ref. ^ and the current work, 
the doping dependence of the 500 meV peak intensity 
is plotted in Fig. [S] For this figure, the peak intensities 
were normalized from sample to sample by the elastic 
(a; — 0) peak intensity. The error bars in relative 
intensity are then large due to the elastic peak intensity 
variations from sample to sample, depending on the 
crystal quality. The general trend, however, is that as 
X increases, the spectral feature intensity decreases, but 
the peak feature survives up to at least a;=0.07, well 
into the superconducting phase. For 2;=0.17, the peak 
appears to be either absent or highly damped. 

It is well known that magnetic correlations are sup- 
pressed upon hole-doping in the cuprates, so that this 
observed doping dependence of the intensity suggests 
that the excitation has a magnetic characteri^ The 
existence of the peak for doping values x = 0.07 shows 
that this excitation survives even if long-range magnetic 
order is absent, so long as significant short-range mag- 
netic correlations are present (which is well known to 
be the case in the superconducting state of LSCO^i^) . 
This is consistent with our observation that the mode 
also survives above the Neel temperature in the undoped 
parent compound. The doping dependence of the RIXS 
intensity is quite similar to the doping dependence of 
the two-magnon Raman spectra in Ref. 37, in which 
there is a strong drop in intensity occurring between 
x=0 and x=0.01. The peak is still visible at x=0.07, but 
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FIG. 6: Peak intensity of the 500 meV peak vs. doping x for 
La2-xSrxCu04. The peak intensities were determined from 
Lorentzian fits of the spectra from Ref. [s^ and in Fig. [O] 
Between different samples, the spectra were normalized by 
the intensity of the elastic line at the same Q position and 
incident energy. The dashed line results from a fit to the form 
I{x) = 7o/(l + ry^), with r ~10. 



no intensity appears on the same scale for x=0.12 (note 
that X in our notation is twice that of Ref. IstI ) . 

Assuming that the 500 meV mode is two-magnon in 
origin, the following physical argument sheds some light 
on the doping dependence of the intensity. The instanta- 
neous spin-spin correlation length ^, which is a measure 
of magnetic correlations, depends on doping as approxi- 
mately ^ ~ a/y/x over a range of dopings, where a is the 
lattice constant.— While there is a certain inverse life- 
time, To, to two-magnon excitations even in the insulat- 
ing compound (energy resolution also contributes to the 
nominal width), we expect the finite magnetic correlation 
length in the doped system to lead to an additional scat- 
tering rate for two-magnon excitations. This scattering 
rate may be calculated by ~ v/^, where v is the char- 
acteristic velocity associated with magnetic excitations. 
Assuming that the spectral weight of the 500 meV mode 
is not a strong function of doping for small x, the peak in- 
tensity is expected to scale as l/(ro-l-ra;) ~ l/{l + ry/x) 
where r ~ However, we caution that using the cor- 
relation length determined in Ref. 70 for low energies 
should not generally apply for these high-energy exci- 
tations. Nevertheless, we find that this functional form, 
I{x) = lo/ i^+ry/x) fits the data quite well up to a; ~ 0.07 
as shown in Fig. [B] Setting Fq — 150 meV consistent 
with the observed Lorentzian linewidth of the mode and 
w ^ Ja as a rough measure of characteristic velocity, we 
find a value for r which is within order of magnitude of 
the value obtained from the best fit to the data, lending 
support to the above argument. 



IV. DISCUSSION 
A. Temperature dependence 

Magnon excitations and magnon-magnon interactions 
in two-dimensional (2D) S=l/2 Heisenberg antiferro- 
magnets have been extensively investigated over many 
years ;^2iiir2^ especially after the discovery of high 
temperature superconductivity. Theoretical calculations 
have shown that the short wavelength magnon near 
the zone boundary is well defined even above the Neel 
temperature^!- that is, the damping due to magnon- 
magnon interactions F(T) is much smaller than the 
magnon energy as long as S> 1. A number of studies 
have looked into the damping of two-magnon excitations 
due to magnon-magnon interactions<^>2&^— The con- 
sensus from these studies is that the magnon-magnon 
interactions alone cannot explain the observed broad 
linewidth of the two-magnon Raman scattering or the 
temperature dependence of the intensity. In order to 
explain the observed linewidth in their two-magnon 
Raman scattering results on EuBa2Cu306, Knoll et. 
al. suggested that magnon-phonon interactionSSiSi is 
an additional source of damping.t55 Such damping from 
phonons is expected to be a larger effect in cuprates 
than in compounds such as K2NiF4, since the magnetic 
energy scale in the cuprates is much larger than typical 
phonon energies, resulting in smaller thermal popula- 
tion of magnons compared to that of phonons i^i^l In 
La2Cu04, the damping may even be more pronounced 
due to a rotational phonon mode which softens towards 
the structural transition temperature at around T=500 

82,83 'pj^jg damping may well be responsible for the 
disappearance of the peak at around the same tempera- 
ture, although further temperature dependence studies 
would be required to test this. 

With the energy resolution used in this study, 
linewidth broadening is difficult to detect. Instead, we 
focus on the shift in the peak energy. Temperature- 
induced softening of the two-magnon energy in 5 = 1 
compounds such as K2NiF4 has been quantitatively ex- 
plained by magnon-magnon interaction theory. Keffer 
and Loudon^i, Bloch®^, and Davies et. al^i^, calcu- 
lated the shifts in the magnon energies nk(r) to be pro- 
portional to the number of magnons determined by the 
Bose distribution. For a Heisenberg antiferromagnet, the 
one-magnon dispersion gets renormalized as 

r!k(r) = a(T)r!k (2) 

q 

where J is the magnetic exchange constant, z is the 
number of nearest neighbors, S is the spin, and the 
rcnormalization factor a{T) needs to be determined 
self-consistently. The constant a(0) incorporates the 
Oguchi®* correction to the spin wave velocity at zero tem- 
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perature. (We note that if the sum in 
nated by a single momentum, it reduce 
Eq. ([T]) which motivated our earUer fit in Q 
sonable to assume that the two-magnoi 
renormaUzed in a similar fashion. In o 
different materials, it is then more usefi 
pression for the energy-shift AE = Eq~ 
magnon excitation, normalized to its 
value -Bo, in the form 

where /(T) is a scaling function obtair [] 

with its argument being the scal( 

f = T/{SJz). Here S is the spin quant 

the number of nearest neighbors, and J 

exchange energy. Such a scaling ansatz i 

interaction corrections to ol{T) are smal 

set a{T) 1 in the Bose function in Eq 

rik cx JzS. To leading order, Eq. Q 

normalized energy shift in terms of the 

parameters of a material. The prefactor of the energy 

shift scales as 1/S since this governs the strength of 

magnon-magnon interactions;-® Such spin-dependence is 

evident from the observed two-magnon Raman energies: 

for spin-1 K2NiF4, the two-magnon energy is renormal- 

ized from 8J to 6.8J, which is a 15% decrease from its 

theoretical value without spin-spin interaction, while 

for spin-1/2 cuprates, the renormalization from 4 J to 

2.7 J represents a 32.5% decrease, or about twice as much. 

We now check if the observed energy shifts can be 
explained by the spin-spin interaction theory by compar- 
ing the shifts in different 2D Heisenberg antiferromagnet 
compounds. The scaling function in Eq. ^ is plotted 
in Fig. [Tja) for our RIXS data in La2Cu04, and the 
two-magnon Raman data in EuBa2Cu30g^ and those 
of K2NiF4.-S2 The peak shifts in the latter material were 
proven to follow the spin-spin interaction theory, at 
least up to r=0.3,®^ so this curve may be considered 
as the theoretical prediction. The clear non-overlap of 
the curves in Fig. El^a) suggests a departure from the 
spin-spin interaction model although the discrepancy 
in the case of EuBa2Cu306 is less than that of La2Cu04. 

There are several reasons why one might expect devia- 
tions from our assumption that all of these materials are 
described by a single-layer Heisenberg model. The pres- 
ence of single-ion anisotropics in K2NiF4, the presence of 
bilayers with strong interlayer couplings in EuBa2Cu306, 
and the effect of temperature dependent phonon-induced 
renormalizations of the antiferromagnetic exchange cou- 
pling, could all plausibly account for the non-overlap of 
the data. Empirically, we find that scaling the tempera- 
ture by the Neel temperature. Tat, appears to provide a 
good coincidence of the peak shift data on all three ma- 
terials as shown in Fig. [Tljb), where T/v replaces JSz as 
the temperature scale. That the data nearly collapse into 
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FIG. 7: (Color Online) Comparison of scaled peak energy 
shift vs. scaled temperature in La2Cu04 as measured by 
RIXS in this study (circles), and two-magnon Raman scat- 
tering in EuBa2Cu306 (triangles) obtained from the peaks 
of the fit curves in Ref. 5^, and K2NiF4 (unfilled squares) 
from Ref. [s^. (a) The peak energy shift is normafized by tiie 
low-temperature energy, and further multiplied by spin S, as 
discussed in the text. The temperature scales are normalized 
by JSz. J of 143 meV is used for La2Cu04, 100 meV for 
EuBa2Cu306^- and 9.7 meV for K2NiF4-^ respectively The 
lines result from fits to Eq. ([T]). (b) Same as (a), except nor- 
malizing the temperature to the Neel temperature Tjv, which 
is 320K for La2Cu04, 430K for EuBa2Cu30a^ and 97K for 
K2NiF4^ (in Ref. [S^ the temperature is already normalized 
to Tn) . 



one curve is strongly suggestive of a spin-spin interaction 
mechanism for the energy shift of the RIXS peak, which 
in turn hints at the magnetic nature of this peak. One 
possible explanation of this data collapse upon using T/v 
scaling is that the 3D Neel temperature already takes 
into account various material specific properties men- 
tioned above which affect the magnetic excitation spec- 
trum. This empirical observation needs to be tested in 
other insulating quantum magnets. 



B. Excitation peak position 

Although the temperature dependence appears to 
confirm that the 500 meV excitation shows characteris- 
tics of magnetic excitations, the question still remains 
as to the relatively weak nature of the magnon-magnon 
interaction within the two-magnon excitation, as well as 
the precise value of the expected q=(7r 0) peak position 
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as determined in units of J, since not all calculations 
yield the same results (see Fig. [31[a)). There are at least 
three possibilities for explaining the observed energy 
in the framework of magnons. Sugai et. al. observed 
a distinct feature in the Raman spectrum of La2Cu04 
at ^550 meV, which was attributed to a 4-magnon 
excitation)^, or alternatively a phonon-induced magnon 
peaki2^i2^ The matter is not resolved,^"' but all scenarios 
for the Raman shoulder-feature involve multi-magnons, 
so much of the temperature and doping dependence 
analysis in this study should still apply. However, the 
dispersion observed with RIXS would suggest that the 
excitation is a two-magnon (unless a 4-magnon has 
similar dispersion, but we cannot find such calculations 
in the literature). We proceed below to explain the high 
energy, with the assumption that the RIXS excitation is 
indeed a two-magnon. 

The 500 meV peak is only observable when the 
incident polarization is along the crystalline c-axis, 
leading to the suggestion in Ref. 'ss' that adjacent ab 
planes could be involved, with the two magnons being 
created on adjacent planes, rather than on a single 
Cu-0 sheet. The small value of J_l would then explain 
the relatively weak interaction between the magnons. 
Donkov and Chubukov also discussed the possibility 
that interplane interactions caused by the c-polarized 
photon could serve to weaken the magnon-magnon 
interaction. Forte et ali^ discussed another mechanism 
that might reduce the amount of energy renornialization 
at the zone-boundary. They found that including 
the effect of longer range couplings (second and third 
nearest neighbors, as well as ring- exchange) increased 
the two-magnon energy at finite q and their calculated 
energy is in agreement with this experiment, as seen in 
Fig. [IJa). Here, we present another possible resolution, 
based on a phase-space argument, which suggests that 
the interaction between in-plane magnons is a strong 
function of the total momentum q. 
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FIG. 8: (Color Online) (a): The two-magnon density of states 
for noninteracting magnons with a total energy u and two 
different total momenta, (b): Intensity plot showing the sin- 
gle magnon states that contribute to give a total momentum 
qtot = (0,0) at uj=ujpeak='^J (where the two-magnon density 
of states is a maximum). The two individual magnons have 
the maximum probability of carrying momenta q2=-qi=k, 
where k lies on the magnetic zone boundary, (c): Same as 
in (b) but for qtot = (7r,0) and ctj=tjpea*:=2 J\/3- In this case, 
the maximum probability is for individual magnons to have 
qi=q2=k, where k=(±7r/2,0), (0,±7r/2). 



We begin by calculating the non-interacting two- 
magnon density of states (DOS). In Fig. [SJ^a), the 
q-integrated DOS is plotted as a function of energy 
with the total momentum qi -f c\2 fixed, where qi and 
(\2 denote individual momenta of the two constituent 
magnons. For the case of qi-l-q2=0, which is relevant 
for two-magnon Raman scattering, the peak in the 
non-interacting two-magnon DOS occurs at an energy 
of 4J, and it is well known that strong magnon-magnon 
interaction brings this down to about 2.75 J. ^^'^^'^^ 
However, the same DOS at qi + q2 = (7r, 0) exhibits a 
broad peak at an energy of around 3.5J (consistent 
with previous calculations, Ref. [s^ for example) with 
a sharp decrease in the DOS above 4J. Even with 
some uncertainty in the value of J, the energy of the 
500 meV peak observed at q=(7r 0) is much closer to 
the non-interacting DOS peak, than is the 2-magnon 
Raman q=(0 0) peak to its corresponding DOS peakfSS, 



This naturally leads to the question of momentum 
dependence of magnon-magnon interactions within a 
2-magnon. 

To gain insight into the weakened interactions in 
the RIXS data, we have looked at the momenta of the 
component magnons, qi and q2 for various total mo- 
menta qi + q2 . Figure IHb) shows constant energy cuts 
of the single magnon contributions to the two-magnon 
density of states at an energy of 4 J and a momentum of 
qi + q2=0. The largest contributions come from single 
magnons carrying the momentum shown as the bright 
color in this figure. One can see that the zone-boundary 
magnons dominate the contributions to the peak in the 
two-magnon density of states. Because the magnon 
group velocity is zero at the zone-boundary, there is 
ample time for these magnons to interact with one 
another. Furthermore, there are a large number of states 
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corresponding to magnons with opposite momenta living 
on the zone boundary into which scattering can occur. 
As a result, interaction effects would be expected to 
be large at this total momentum - as is observed in 
the Raman experiments. In contrast, if we look at the 
single magnon contributions at a total momentum of 
qi + q2 = (7r, 0) and an energy of 3.5 J, Fig. IDJc), we find 
that the individual magnons contributing to this peak 
are not concentrated in the zone boundary region. For 
example, a pair of magnons near the (7r/2, 0) position 
could add up to form a two-magnon mode with total 
momentum (tt, 0). Since these magnons can generally 
have a nonzero velocity with respect to each other we 
argue that the effective interaction between these two 
magnons is diminished, providing a natural qualitative 
explanation for the observed excitation energy. We 
note that a similar idea was also proposed by Vernay 
et. al., who argued that magnon-magnon interaction 
is weakened because the average distance between 
spin-flips is larger when the total momentum is finite. ^« 



V. CONCLUSION 

We report a comprehensive investigation of the 500 
meV peak in La2-xSrxCu04 observed with Cu i^T-edge 
resonant inelastic x-ray scattering (RIXS). We have 
carried out studies of incident energy, momentum, 
temperature, and doping dependence. At temperatures 
below 300 K the peak energy shift of the RIXS peak 
is quantitatively consistent with both Raman-'^^- and 
neutron^'' studies of La2Cu04. At high temperature 
the precipitous drop in the peak intensity is suggestive 
of magnon-phonon interaction. However, the softening 



of the mode appears to be a result of coupling to 
spin fluctuations, as indicated by the scaling with T^r. 
These quantitative and qualitative similarities with two- 
magnon Raman scattering, together with the observed 
dispersion and doping dependence, which are consistent 
with the expectations for a two-magnon excitation, 
provide compelling evidence for the magnetic nature of 
the 500 meV peak in La2-xSra;Cu04. 
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Estimating effective magnetic exchange constant J (includ- 
ing the Oguchi correction) based on the 2-magnon Raman 
peak calculation of 2.75J would imply that the 500 meV 
peak is 3.44J. On the other hand, a neutron scattering 
study by Hayden et al. [Hayden et al., Phys. Rev. Lett. 67 
3622 (1991)] yielded an uncorrected exchange value of 132 
meV, which was at room temperature. Assuming a 4 per- 
cent increase of J upon going to low temperaturCj^'Z the 500 



meV peak energy would then be 3.12J (using a spin-wave 
velociy correction factor of 1.17), which is 11 percent be- 
low the unrenormalized energy of 3.5J for q=(7r 0). While 
this would represent a significant decrease, it is still con- 
siderably less than for the 2-magnon Raman q=(0,0). For 
the same value of J, the 2-magnon Raman energy would 
be 2.5J, which is a 37 percent decrease from its nominal 
value of 4 J. 

®® When they are exactly at (7r/2 0), the relative velocity is 
zero, but the weight of such cases would be small. 



Appendix A: Background Subtraction Methods 

Due to thermal and static diffuse scattering, there 
exists a large quasi-elastic (a; = 0) contribution to the 
background, and the tail of this elastic background 
makes it difficult to study excitations below ^500 
meV. For this reason, it is important to subtract 
this quasi-elastic background to obtain the excitation 
spectrum. We have employed three different methods to 
determine the background, depending on experimental 
convenience or circumstance. The first is to utilize the 
high-temperature spectrum as background. This is most 
effective for studying temperature dependence, and is 
commonly used in many other spectroscopy experiments 
and is described in detail in Sec. IIII Bl The other two 
methods are specific to RIXS, and it is worth discussing 
these in detail here. 

In the first of these, the "off-resonance" spectrum - 
measured with an incident photon energy well below the 
resonant energy - is subtracted from the "on-resonance" 
spectrum, obtained with the incident energy set to 
slightly above the absorption peak. The off-resonance 
spectrum is normalized such that the intensity of the 
quasi-elastic tail on the energy-gain side best matches 
with that of the on-resonant spectrum, as shown in 
Fig. [n{a)-(b). It is necessary to do so because the 
absorption changes dramatically in this energy range, 
resulting in a large change of the quasi-elastic intensity. 
While this method utilizes the unique resonance prop- 
erty of RIXS excitations, and eliminates artifacts from 
asymmetric analyzer tails, it can be time consuming 
to obtain both on- and off-resonance scans with good 
statistics. 

The final method relies solely on detailed balance. 
Since we are measuring excitations at a much higher 
energy than the thermal energy, the intensity of an 
excitation on the energy gain side (w < 0) is essentially 
zero due to the principle of detailed balance. Therefore, 
the energy gain spectra represents contributions purely 
from the resolution function and the background. Thus, 
subtracting a symmeterized version of the energy gain 
spectrum should remove both contributions, leaving just 
the inelastic processes. This is accomplished by reflecting 
the spectrum about a; = (cj — > — w), and subtracting 
this energy gain spectrum from the original energy loss 
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FIG. 9: (Color Online) (a) The RIXS intensity versus pho- 
ton energy loss, for Lai.93Sro.o7Cu04 at q=(7r,0) and T=20 
K, for incident energies "on-resonance" 15^=8993 eV (filled 
squares) and "off-resonance" (empty circles). For improved 
statistics in the "off-resonance" spectrum we combined spec- 
tra of i5i=8981 eV with i5i=8987 eV, whose low-energy spec- 
tra were identical within error bar. The resultant spectrum 
was then normalized to match the energy-gain side of the "on- 
resonance" spectrum, (b) Subtraction of the spectra shown 
in (a), which has reasonable error bar only above 0.3 eV en- 
ergy loss, (c) The same "on-resonance" spectrum plotted as a 
function of energy loss (filled squares) and energy gain (open 
circles) , and a fit of the energy gain spectrum to a Lorentzian 
lineshape plus constant (filled diamonds), (d) Subtraction of 
a fit of the energy gain spectrum from the original spectrum. 



spectrum for w > 0, as shown in Fig. Etc)-(d). This 
method assumes that the resolution function (dominated 
by the analyzer) is symmetric. 

Due to the steep tail of the elastic peak, the subtrac- 
tion of the background spectrum by these methods can 
be quite sensitive to even small energy shifts 5u of ana- 



lyzer. Small temperature variations in the experimental 
hutch, for example, cause the lattice parameter of the 
analyzer to change, resulting in a 55 meV/°C shift of the 
actual energy. The largest variations, of up to a couple 
of degrees °C, can typically occur after the hutch is first 
closed. To compensate for this effect, the temperature 
near the analyzer was recorded for each data point and 
the energy values adjusted accordingly during the 9ID 
measurements. Apart from the first several hours after 
the hutch is closed, the typical drift over the course of a 
scan was of the order of 20 meV or less, which mostly 
affected the part of the spectrum below 0.4 eV. 

While these methods of estimating the contribution 
from the elastic tail should be equivalent in ideal 
conditions, a comparison of Fig. [9] (b) and (d) reveals 
that the energy-gain side subtraction method can result 
in an additional sloping background when compared 
to the off-resonance method. (It is important to note, 
however, that both methods give the same value for 
the peak position, the parameter of importance in this 
work.) This might be due to non- resonant inelastic 
scattering which would also be present in (and thus 
canceled by) the off-resonance scan. Such non-resonant 
scattering could be due to intraband excitations, and 
are thus expected to be more prevalent in doped samples. 

In the course of these experiments, different circum- 
stances favored one method or the other (for example, 
if the elastic peak happened not to be perfectly sym- 
metrical about uj=Q due to analyzer imperfections, 
this would be problematic for the energy-gain side 
subtraction method). For the 9ID and BLllXU data 
we used off-resonance scans for background, while the 
energy-gain side subtraction method was employed in 
the 30ID momentum dependence scans. It should be 
emphasized that the results from using the different 
background subtraction methods were found to be 
consistent with each other and generally fall on the same 
curves for temperature and momentum dependence (the 
one exception is discussed in llll A[) . 



